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Abstract 

The tendency of amphiphilic molecules to form micelles in aqueous solution is a consequence 
of the hydrophobic effect. The fundamental difference between micelle assembly and macroscopic 
phase separation is the stoichiometric constraint that frustrates the demixing of polar and hy- 
drophobic groups. We present a theory for micelle assembly that combines the account of this 
constraint with a description of the hydrophobic driving force. The latter arises from the length 
scale dependence of aqueous solvation. The theoretical predictions for temperature dependence 
and surfactant chain length dependence of critical micelle concentrations for nonionic surfactants 
agree favorably with experiment. 
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I. INTRODUCTION 



This paper concerns the formation of micelles, which are the simplest form of amphiphilic 
assemblies. Our treatment of this phenomenon is based upon the length scale dependence of 
hydrophobic effects . Namely, the free energy to solvate small hydrophobic molecules 
scales linearly with solute volume, while that to solvate large hydrophobic species scales 
linearly with surface area. The crossover from one regime to the other occurs when the oily 
species presents a surface in water extending over about lnm 2 . Due to these contrasting 
scalings, the free energy to solvate a collection of small oily species that are well separated 
in water can exceed that of solvating a large cluster of these same species. The resulting free 
energy difference is the hydrophobic driving force for assembly. While this force that derives 
from this scaling difference is surely the mechanism that drives oil-water phase separation, 
one may wonder if it is also applicable in the case of amphiphilic assembly where interfaces 
between solute clusters and water contain polar or charged head-groups as well as oil. In 
this paper, we argue that it is indeed applicable. 

The primary difference between oil-water phase separation and amphiphilic assembly is 
due to stoichiometry. Each amphiphilic molecule contains an oily species that is constrained 
to remain within a molecular length of a hydrophilic species. For large clusters of am- 
phiphiles, this constraint leads to an entropic cost for clustering oily components that grows 
faster than cluster volume. Thus, unlike simple oil-water phase separation which forms 
macroscopic domains, growth of amphiphilic assemblies is limited to mesoscopic domains. 
Our treatment of this effect uses an estimate of entropy ^ obtained from an electrostatic 
analogy for stoichiometric constraints jj, 0, |(| . In the next section, we describe this estimate 
along with the other factors that contribute to micelle formation. The resulting expression 
for the critical micelle concentration, p cmc , is then compared with experiments in Section 
III. In particular, we show that our expression yields good predictions of this concentration, 
changing with temperature and with amphiphile chain length in accord with experimental 
observations. An appendix is used to augment the discussion in Section II. 
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II. THEORY 



A. Law of mass action 

We consider an aqueous solution of neutral amphiphilic molecules (i.e., non-ionic surfac- 
tants), each of which has a single alkyl chain as its hydrophobic tail. In general, amphiphiles 
can form aggregates of various sizes and shapes. We will assume each micelle is spherical 
and neglect the effects of fluctuations in micelle size and shape. Thus, we imagine that each 
surfactant molecule exists either as a monomer or as part of a spherical n-mer. We denote 
the number densities of the monomers and n-mers by pi and p n , respectively, so that the 
total surfactant concentration is given by p = p\ + np n . 

The concentrations of monomers and micelles are related by the law of mass action 

p n a 3 = ( Pl a 3 ) n exp(-/3AG), (1) 

where j3 denotes inverse temperature (i.e., f3~ l = k^T), a is a microscopic length that 
specifies the standard state convention, and 

AG = f n -nf 1 (2) 

is the driving force for assembly, namely, the free energy of the n-mer, /„, relative to that 
of n monomers, nf\. We take a to be approximately the girth of a surfactant molecule, see 
Fig. El 

For large n, eq^hnphes the existence of a threshold concentration of surfactant molecules 
p cmc , at which the density of aggregates becomes significant. Because this crossover is 
precipitous, its location is almost independent of the specific definition of the threshold as 
long as it is physically sensible. Specifically, to within corrections of order n~ l Inn, 

\np cmc a 3 = f3AG/n*. (3) 

The driving force per surfactant, AG/n, is a function of n, and it is to be evaluated at the 
most probable aggregation number, n*. This number is the value of n that minimizes AG/n. 
Equation |3] is discussed further in the Appendix. 
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B. Driving force 



The contributions to AG can be found in three steps, employing the thermodynamic 
cycle illustrated in Fig. [TJ 

1) Creation of a cavity. A micelle will fill a region vacated by water. Assuming the extent 
of the surface is at least 1 nm 2 , the free energy to create this cavity is 

AG 1 = 7 A (4) 

where A denotes the surface area of the cavity, and 7 is the water-vapor surface tension. 
In general, there is also pressure-volume work for forming a cavity in a liquid. For water 
at standard conditions, pressure is sufficiently small that this contribution is negligible for 
cavities with diameters less than 5 nm. We will limit our consideration to sizes within this 
range. 

2) Filling the hydrophobic core. Imagine disconnecting each hydrophobic tail in a surfac- 
tant from its respective hydrophilic head group and moving the hydrophobic tail from water 
into the micelle core (we will reconnect the heads and tails in Step 3). A total of n tails 
must be moved to fill the cavity formed in Step 1. As such, one part of the free energy to fill 
the cavity is — n A/i, where — A/z is the free energy change in transferring the hydrophobic 
tail (e.g., an alkane chain) from water into the oily hydrophobic core. An additional part of 
the free energy for filling the cavity is an interfacial contribution due to the presence of van 
der Waals attractions between oil and water. These interactions cause the oil-water surface 
tension, 7 ow , to be lower than the water- vapor surface tension, 7 (see, for instance, Ref. 
Thus, the free energy for filling the cavity is 

AG 2 = -n An - A 7 A, (5) 

where A7 = 7 - 7ow . 

The interior of a micelle is densely packed and much like a hydrocarbon liquid |i| . Thus, 
A/j, is close to the transfer free energy for moving the associated alkane chain from oil 
into water. However, it is slightly smaller than this value because the environment of an 
alkane chain in a micelle interior is more confining than that in bulk oil [9]. The numerical 
consequence of this small difference will be discussed in Sec. III. To the extent that the 
micelle is spherical, A = 47rL 2 , where L is the micelle radius. Since the interior is densely 
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packed, L is given by 47rL 3 /3 = nSa 2 , where S is the mean length over which a polar head 
group is separated from an alkyl group within a surfactant molecule, see Fig. 1. From these 
considerations, 

AG 1 + AG 2 = -nA^i + gn 2 ^, (6) 

where g = (36tt) 1 / 3 ( 7ow a 2 )(5/a) 2 / 3 « 4.8 x ( 7ow a 2 )(5 / a) 2 / 3 . 

Th e rig ht-hand side of eq El is essentially the free energy for nucleating oil clusters in 
water jlOj]. It is the hydrophobic driving force identified in the Lum- Chandler- Weeks the- 
ory Q| . The first term is proportional to the volume of hydrophobic units. The second term 
is proportional to the area of the interface. The first term is extensive in n and dominates 
at large n. Thus, if only AG\ and AG 2 were significant, the strength of the driving force 
would grow without bound leading to macroscopic clusters. However, these contributions 
are counter balanced by a third term that we consider now. 

3) Placing hydrophilic head groups on micelle surface. In the final step, the hydrophilic 
head groups are reconnected to the hydrophobic tails, placing them at the water-oil interface 
so as to maintain a favorable solvation energy. This positioning is to be done while simulta- 
neously enforcing the connectivity between heads and tails and while also maintaining the 
densely packed interior. These conditions result in an entropic cost that increases super- 
extensively with aggregate size. The form of this third contribution to the driving force 
is conveniently estimated from the electrostatic analogy of stoichiometric constraints jj, 0] . 
The result is jsj] 

AG 3 = hn 5/3 /P, (7) 

where h = (3/(47r)) 2 / 3 (96/49)(a/<5) 4//3 ps 0.75 x (a/5) 4//3 . In employing this analogy, it is 
important to note that the micelle volume is essentially that of the densely packed alkyl 
chains. 

C. Micelle size and critical micelle concentration 

Combining the three contributions discussed above gives the driving force in units of k^T: 

j3AG ^ -n(3 Afi + (3 g n 2/3 + h n 5/3 . (8) 
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Minimization of AG/n therefore gives 

n* « /3g/2h = (49tt/48) /3 7 <5 2 . (9) 
With this aggregation number, eqs El and El yield 

In p cmc a 3 = c (/?7 ow a 2 ) 2/3 - /?A/x, (10) 
where c = (5832/49) 1 / 3 ss 4.9. 



III. COMPARISON WITH EXPERIMENT 

Equation El is the principal result of this paper. It expresses the critical micelle con- 
centration in terms of measured quantities and one adjustable parameter, the molecular 
length scale a. As such, its validity is easily checked. Here we do so by considering m- 
alkyl hexaoxyethylene glycol monoethers, C m E 6 , i.e., (CH 3 )(CH2) m _i(OCH 2 CH2)60H. The 
critical micelle concentrations for this class of nonionic surfactants have been determined 
experimentally. 

The surface tension is required to compare eq El with experimental results. At room 
temperature, the oil-water surface tension is 7 ow ~ 51mN/m Were we to neglect 

A7, we would instead take the room temperature value of the water-vapor surface tension, 
7 = 72mN/m The difference A7 = 7 — 7 ow is mainly enthalpic As such, a good 
approximation for the temperature dependence of the required surface tension is d^ OVf /dT w 



d-f/dT « -0.17mN/m/K 12]. 



The transfer free energy, A/x, is also needed. To the extent that the micelle interior is 



like a bulk hydrocarbon liquid, it can be obtained from solubility measurements 13, 14]. 
This free energy change for transferring alkane chains with fewer than 12 carbons from oil 
to water depends linearly on the number of carbons, m. The linear fit to experimental data 
is A/i (m) = (2.25 + 0.9 m) kcal/mol at room temperature. We have used the subscript "0" 
to indicate that this free energy for transferring between bulk phases must differ to some 
extent from that for transferring from a micelle interior to water. Indeed, experimental 
evidence for alkane chains with m < 6 9] indicates that A/x(m) w (1.9 + 0.77m) kcal/mol = 



A/x (m) + AA/x(m). 

The temperature dependence of A/x can be estimated from the observation o 



temperatures for transfer entropies and enthalpies, Tg and T#, respectively 
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convergence 
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Namely, 



A/,(T,m) « AC p (m) [(T - T*) -T\n(T/T*)} , (11) 

where T£ = 112 C and T# = 22 C [17]. The heat capacity AC P (m) is chosen so as to fit the 
experimental A/i(m) at room temperature, i,e., 

AC p (m) = [A^ (m) + AA/j(m)] / [(298K - T£) - 298K ln(298K/T^)] . (12) 

EquationEltogether with the dependence of transfer free energies on alkyl chain predicts 
that p cmc varies exponentially with m. When plotted on a logarithmic scale, it should have 
the slope —dAfi/dm. Figure El compares this prediction with experimental data. The 
measured data is indeed linear over a range of concentrations spanning several orders of 
magnitude. The slope of the experimental data is approximately -0.7 kcal/mol, which is 
close to the value of -0.9 kcal/mol inferred by neglecting AA/i(m), and even closer to the 
value of -0.77 kcal/mol inferred by accounting for this correction to transfer free energies 
between entirely bulk phases. 

Figure E3 concerns the temperature dependence of the critical micelle concentration. The 
experimental data for C^Eg show a minimum in the critical micelle concentration near 50 
Celsius. A similar minimum is found from eqEl The excellent fit to experimental data is 
obtained by choosing a so as to match the experimental data at T = 25 °C. This method 
of choice gives a = 3A, which is not an unreasonable value for the microscopic length. 
Neglecting the corrections A7 and AA/i lead to similar lengths, and similar though inferior 
fits to the data over the observed temperature range. 

Since f3^a 2 is of order 1, and 8/ a is of order 10, we see that this expression predicts 
aggregation numbers of order 100. This cluster size is consistent with our assumption of 
large n. 



IV. DISCUSSION 



Elements of the theory we have presented here for micelle assembly can be found in 
earlier works. For instance, the surface energy and entropy terms, AGi and AG3, have 

As a result, that paper obtains the same scalingof aggregate 



been considered in Ref. 



size as given by eqlHl namely, n* oc f3^5 2 . By neglecting AG2, however, Ref. errors in its 
treatment of the temperature dependence of p cmc . 



Detailed consideration of the transfer free energy contribution AG2 and its linearity 

n 

with m is found in Tanford's monograph [9j]. That work does not consider AGi and AG3 
which compete with AG2. Both AGi and AG 2 contribute significantly to the temperature 
dependence of p cmc . 

For the neutral surfactants we have considered, the net hydrophobic driving force, 
AG 1 + AG2, is balanced by the super-extensive entropy term, AG3. It is remarkable how 
well p cmc is described over a wide range of conditions with only these three terms. For 
ionic surfactants, an additional super-extensive term will limit the micelle size. It is an 
electrostatic contribution that opposes the clustering of like charge. For low concentrations 
of counterions, its scaling with n is identical to the entropic term, growing as n 5//3 for the 
case of spherical micelles. The size of this term will depend upon ionic strength in a fashion 
not yet determined. 
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APPENDIX 

Here we consider further the law of mass action, considering the water-surfactant mixture 
a 8 an idea, solution of Octant, and Octant „tes. The fr ee ene rgy density is Q 

PF = ^ [Pn (hip n a 3 - l) + pnPfn ~ PPS^pn] ■ (13) 
n 

The first term accounts for the translational entropies of the aggregates. f n is the internal 
free energy of an n-mer, and /is is the surfactant chemical potential that controls the total 
surfactant density p = J2 n n Pn- 

The equilibrium partitioning of densities is that which minimizes eq ED 

p n a 3 = exp {-(3f n + n(3ps) (14) 

from which the law of mass action, eq^ follows. 
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To the extent that micelles are mono disperse, p n ~ P i5 n ,i + Pn*b~n,n*- Here, 5 ntn * refers to 
the Kronecker delta. Substituting this expression into eq IT31 gives 

(3F^ Pl (ln Pl a 3 -l)4 P ~ 



n 



In I ^a 3 



n 



1 



+ PPlfl + P [ P —? L ) /n—/W (15) 



In this approximation, the equilibrium partitioning of aggregated and unaggregated surfac- 
tants is obtained by minimizing this expression with respect to p\ and n* . That minimization 
gives again eqUJ and also 



= df3F/dn* 



d/3AG/n 



- (lnpia 3 - f3AG/n) 



(16) 



dn n 

The critical micelle concentration is identified as the lowest surfactant density at which a 
measurable fraction of surfactants, x, aggregate to form micelles. As such, np n = xpi, and 
the law of mass action is exp (— (3 AG) = {xpia 3 /n) /(pia 3 ) n . Therefore, 



pAG/n = (1 - 1/n) \n Pl a 3 + (1/n) In (n/x) 



(17) 



and p cmc = pi [1 + O (x)]. Equation El with the condition that n* minimizes AG/n thus 
follows from eqs ITBI and IT71 when n is large. 

An alternative definition of the critical micelle concentration is the lowest density oi 
at which eq El has a non-trivial solution for n*. This convention is used in Ref. |3|. 
The free energy develops a local minimum for n > 1 at the concentration pa 3 = 
exp d (/?7o W o 2 ) 2/3 - (3Ap , where c' = (4320/49) 1/3 rj 4.45, and the location of the min- 
imum is at n* = j3g/ (5h) . Comparison with eqs El and El respectively, shows that this 
convention is in close accord with the approach we have taken. 
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FIG. 1: Thermodynamic cycle of micelle formation: the process of assembling n separated am- 
phiphiles (a) to a micelle (d) can be performed in three steps: 1 ) Creating a cavity in the solvent 
(light gray) (b); 2) Transferring the hydrophobic chains (dark gray) from the aqueous solution into 
the cavity (c); 3) Distributing the polar units (gray) over the surface of the cavity, and reconnecting 
them to the hydrophobic groups (d). 
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FIG. 2: Critical micelle concentration of C m Eg as a function of the length of the 



rydrophobic tail 

at room temperature. Crosses and circle are experimental data from Refs. iTjJ and 12(1 respectively. 
A reference length ao ~ 3 A was chosen to transform the reported values into dimensionless units. 
The curves show the results of eqEIl using the transfer free energy Afi (solid line) and the bulk 
approximation A/io (dashed line). The length a is chosen such that the experimental value at 
m = 12 from |2(| is recovered. Solid line: a/ao = 1, dashed line: a/a® = 1.3, or a/ao = 1.6 if the 
correction A7 to the surface tension is neglected. 
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FIG. 3: Temperature dependence of the critical micelle concentration of C12E6. The meanings of 
the symbols are the same as in Fig. |2I The solid line is the prediction of eq 1101 Disregarding the 
correction AA^ yields the dotted line, and additionally neglecting A7 results in the dashed line. 
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